Mon. Not. R. Astron. Soc. 000, [T]-?? (2007) Printed 3 February 2008 (MN UT^ style file v2.2) 



Host Galaxy Extinction of SNIa: Co-evolution of ISM 
Structure and Extinction Law with Star-Formation 

B. W. Holwerda^* 

''^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 



Accepted ?. Received 29 dec 2007; in original form 



ABSTRACT 

This paper presents a mechanism that may modify the extinction law for SNIa ob- 
served at higher redshift. Starting from the observations that (1) SNIa occur predom- 
inantly in spiral galaxies, (2) star-formation ejects ISM out of the plane of spirals, (3) 
star-formation alters the extinction properties of the dust in the ISM, and (4) there 
is substantially more star-formation at higher redshift, I propose that spiral galaxies 
have a dustier halo in the past than they do now. The ejected material's lower value 
of Rv will lead to a lower average value (Rv) for SNIa observed at higher redshift. 

Two relations in SNIa observations indicate evolution of the average Rv- the re- 
lation of observed Ry with inclination of the host galaxy at low redshift and the 
matching of the distribution of extinction values {Ay) for SNIa in different redshift 
intervals. The inclination effect does point to a halo with lower Ry values. In con- 
trast, the distributions of Ay values match best for a Rv{z) evolution that mimics 
the relation of SNIa dimming with redshift attributed to the cosmological constant. 
However, even in the worse case scenario, the evolution Ry can not fully explain the 
dimming of SNIa: host galaxy extinction law evolution is not a viable alternative to 
account for the dimming of SNIa. 

Future observations of SNIa -multi-color lightcurves and spectra- will solve sep- 
arately for values of Ay and Ry for each SNIa . Solving for evolution of Ry (and 
Ay) with redshift will be important for the coming generation of cosmological SNIa 
measurements and has the bonus science of insight into the distribution of dust-rich 
ISM in the host galaxies in the distant past. 

Key words: (cosmology:) distance scale cosmology: observations galaxies: high- 
redshift galaxies: ISM (ISM:) dust, extinction (stars:) supernovae: general 



uncertainty in SNIa observations because the applicable ex- 
tinction law remains poorly understood. This will need to be 
addressed in order to use SNIa in the next step in accuracy 
for a cosmological probe. Dust attenuation affects both the 
observed SNIa rate and the distance modulus. Extinction 
by dust occurs in three instances before observation of the 
SNIa light: (1) in our own Milky way, (2) in intergalactic 
space, and (3) in the host galaxy of the SNIa. 

Galactic extinction (1) i s a well-studied p r oblem , be- 
cause it is an ubiquitous one. iBurstein fc HeilesI l|l984 ) pro- 
duced a map of Galactic extinction based on HI maps 
an d distant galaxy cou nts. It was superseded by the map 
of ISchlegel et al.1 l|l998h . based on CORE and IRAS far- 
infrared and sub- mm maps. The latter map is now gener- 
ally used, together with th e Galactic Extinction Law to cor- 
rect extragalactic sources l|Cardelli et al.l Il989l : iFitzpatrickl 
I1999D : the inferred extinction along a line-of-sight is pro- 
portional to the reddening: Ay — Rv x E(B — V). The 
canonical value of Ry is 3.1 with occasionally lower values 
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1 INTRODUCTION 

Supernova type la (SNIa) distance modulus measurements 
have grown into a powerful measurement of the equation 
of state of our Universe. Accurate cosmological distances 
combined with redshift measurements allow for a precise 
characterisation of the Hubble flow as well as the addi- 
tional acceleration attributed to the cosm ological constant 
jRiess et al.lll998l : iPerlmutter et al.1 ligggl ). The increasing 
statistics of SNIa measurements together with more infor- 
mation for each separate supernova event have progres- 
sively lowered the observational uncertainties: dust extinc- 
tion , photometric eirorM 



e.g.jTonrv et al.ll2003l: iKnop et al.ll2003l : iBarris et al.ll20o4 
IConlev et al.ll2006l : lAstier et al.ll2006D . 

Extinction by dust remains a problematic systematic 
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toward s the Galactic Center (lUdalski 20031) a nd higher else- 
where l|Cardelh et al]ll989l : iFTtzpatrickHlOOgl ). 

Extinction by dust in intergalactic space (2) has been 
proposed as an alternative explanation for SNIa dimming, 
which is generally at tri buted to the cosmologica l acceler- 
ation (|Aguirrd [T999al lbl: lAguirre fc Haim"anll2000l ). The re- 
sulting extinction law of this dust would effectively be grejQ 
because the attenuating dust would be spread over all red- 
shifts along every line-of-sight. The coincidence of both uni- 
form distribution in both sky and redshift space does give 
this explanation a somewhat contrived appearance. The in- 
jection of dust into the intergalactic medium would have 
to be constant and substantial. Models of a d usty uni- 
verse (|Goobar et al.ll2002l : [Robaina fc Cepall2007l ) find this 
grey dust explanation increasingly inconsistent with obser- 
vational data (See also e.g., .Riess et al. 2007). 

Extinction within the SNIa's host galaxy (3), dust in the 
immediate surroundings and any disk, ring or spiral arm the 
line-of-sight passes through in projection, is an observation- 
ally evident, yet not fully constrained uncertainty. The Dark 
Energy Task Group Report ( A lbrecht et al. 20061 notes this 
as a primary source of uncertainty for SNIa measurements. 

Three characteristics of the host galaxy's dust could - 
and are expected to- change over the history of the universe: 
(a) total dust content or mass, (b) dust distribution within 
the host galaxy and (c) dust composition. The overall effect 
on the effective extinction law is of interest for the distance 
determination from SNIa light curves. 

Dust mass (a) is a variable in several Spectral 
Energy Distribution (SED ) studies of distant galaxies. 
ICalzetti fc HeckmanI l|l999l ) modeled the overall dust con- 
tent of galaxies over time and fo und that a maxim u m oc - 
curred either at z=l or at z=3. 'Rowan- RobinsonI (120031 ) 
modeled the SED's from distant galax ies and similarly found 
a ma ximum dust content at z=l. Ilglesias-Paramo et al.l 
l|2007t ) find a steady increase of dust mass with time from 
the UV-IR SED of gala:xies. The typical dust mass found 
in distant galaxies is very much a function of the se- 
lected sample. Far-infrared selected samples point to dust- 
rich galaxies, similar to Arp 220 (|Rowan-Robinson et al.l 
I2OO5I ). optical/UV select ed samples point to disks very sim- 
ilar to the local ones (jSaiina et al.l I2OO6I) and Lyman- g 
galaxies point to low-extinction disks ( Nilsson et al.|[2"007^ . 
However, more dust mass should not affect the SNIa dis- 
tance determinations if the extinction law remains the same 
for nearby and distant galaxies. More dust in the distant 
galaxies will p redom inantly affect the observed SNIa rate 
(Hatano ct al. 1998; CappcUar o et al. 1 19991: iGoobar et al.l 
[2002,; Riello fc Patat 2005 ; Mannucci et al.ll2007^ ■ as heavilv 



obscured SNIa drop from observed samples. 

The dust distribution in host galaxies (b) is commonly 
modeled as a double exponential, one radial and one verti- 
cal, sometimes with rings to mimic spiral arms. The radial 
scale of the dust distribution is assumed to be similar to 
the stellar one and the vertical dust scale is supposed to 
be much smaller than the stellar one. Previous observations 
of the scale-height in nearby (edge-on) galaxies appeared to 
corroborate the small scale-height fe.g. IXilouris et al]|l999l : 



^ The term grey extinction is used to denote that there is no 
relation between the reddening of an object and its extinction. 



iHolwerda eraLll2005l : lBianchill2007l ) but recent observations 
indicate a much higher scale -height for the dust (|Seth et al.l 
I2OO5I : iKamphuis etal]|2007l ). similar to the stellar scale. If 
the average scale-height of the dust distribution was higher 
in the past, then SNIa in the plane of the disk will encounter 
more extinction, especially when viewed in an inclined sys- 
tem. The different distribution has a similar effect as varia- 
tion the dust content of galaxies and only the observed SNIa 
rate will be affected, unless dust composition and distribu- 
tion are related. 

The dust composition (c), notably the ratio of small 
to large grains, dire ctly affects th e observed extinction law 
(See the review bv lDraind lioO^. Evolution in the aver- 
age extinction law is the most troubling possibility, as this 
would affect the extinction correction of SNIa and indirectly 
the measured Hubble flow and acceleration. For local SNIa, 
variations in the extinction law h ave been obs erved (e.g., 
iRiess et al.lll996l : Ijha et al.|[2007l ). IWand (I2005D explained 
the different observed extinction law for some SNIa as the 
effect of circumstellar material around the SNIa progenitor. 
This is a plausible scenario as the massive-star supernova e 
are very efficient dust producers llSugerman et al] I2OO6I ) . 
iPatat et al.1 (|2007l ) and Iwang fc Mukheried l|2007l ) report 



observations of such material. 

Alternatively, there is substantial evidence for a link 
between star-formation a nd extinction l aw in star-forming 
galaxies (See the review in ICalzetti|[200ll , and the references 
therein.). Star- formation produces more small grains and the 
intense UV fields alter grain composition. 

Since it is not unreasonable to suppose that star- 
formation affects both the distribution and the extinction 
characteristic of dust in spiral disks, I propose a simple 
model for the evolution of the extinction law applicable to 
SNIa measurements. The aim of this paper is to explore this 
link between the extinction law for SNIa and star-formation 
of host galaxies and to investigate how much extinction law 
evolution could reasonably influence SNIa distance measure- 
ments. The paper is organized in a brief review of current 
observations of evolution in host galaxy dust and Rv (§2), 
a description of the model (§3), two tests based on current 
SNIa observations (§4), and discussion and conclusions (§5). 



2 OBSERVATIONAL STATUS 

Observational evidence for a different extinction law at 
higher redshifts comes from lensed quasars, QSO redden- 
ing by damped Lyman-Q systems and GRB afterglows, as 
well as SNIa measurements. 

Gravitational lenses of QSOs find a range of values for 
R y up to a re dshift of 1 (Nadeau et al."l99l'; 'Falc o et al ' 
Motta et al. 2002; Toft et al. 2000; Goicocchea et al 



199c 



20051 : lElfasdottir et al.l |2006) . However, the lenses are of- 



ten ellipt i cal g alaxies, which are not typical SNIa hosts. 
I Pel et aU (Il991 ) find evidence of Q SO reddening by Lyman - 
Q systems but this is disputed bv iMurphv fc Li ske (2003), 
both based on the SDSS sample. lYork et al.l (^2006) find 
a SMC-type exti nction based on a ll SDSS QSOs with a 
Lyman-a system. iKann et al.1 (|2006l ) compare extinciton for 
different GRB afterglows and find evidence for dust in the 
GRB hosts. None of these observations effectively constrain 
the extinction law applicable to SNIa. 
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l2005h . (II) star- formation ejects 
of spiral disks and the ejected 
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Figure 1. A cartoon of our model: star-formation ejects plumes 
of ISM with lower Ry values out of the disk of the host galaxy. 
A line-of-sight to a SNIa can intersect such plumes. The higher 
star- formation at z~l results in more plumes and hence a higher 
probability for a SNIa line-of-sight intersecting one. 



In the SNIa literature, the proble m of SNIa host galaxy 
dust extinction is well recognised jAlbrecht et al. 20061: 
Conlev et ahl l2006l : lAstier et ah] l2006l : IWood-Vasev et ahl 
20071 ). Observations of high redshift SNIa tell us that 



their spec tra are similar to local ones (jGaravini et al.ll2005l : 
iHook et a l. 2005j^ Garavini ^_ct al. 200j), with no anomalous 
reddening l|Knop et al.ll2003i ) and dus t masses of host galax- 
ies si milar to those of local galaxies (IClements et al.ll2004l . 
l2005l . based on sub-mm data). 

Several studies look at the relations between SNIa 
lightcurve properties (peak brightness, duration, and colour) 
and host galaxy properties (type, stellar mass, extinction, 
and star- formation). A relation between SN peak bright- 
ness a nd ho st galaxy type is we l l estab l ished (plamuv et ahl 
1996', 2000'; Saha ct ajj Il997l. 'l999: Parodi ct al.' '200d!; 



Sandage et al 2001: Sullivan et al. 2003 : Reindj et a l. 2005i). 
Sullivan et ahl ( 20061 ) refine the relation to one between the 



SNIa lightcurve and star-formation in the galaxy: passive 
galaxies, with no star-formation, preferentially host faster- 
declining/dimmer SN la, while brighter events are found 
in systems with ongoing star- formation. The SNIa related 
to ongoing star-format ion could well dominate the popula- 
tion at higher redsh i ft |Sullivan et al.ll2006l : iMannucci et ahl 
l2007l : iHowell et al. I l2007l ). as there is much more star- 
formation. Host galaxy extinction is of equal importance 
as any change in the SNIa population. ISulhvan et ahl (|2003l ) 
report dimmer SNIa for spiral galaxy hosts, especi ally when 
rejec ted in or through the inner part of the disk. lArbutinal 
2007t l also finds a radial dependence of the extinction of lo- 
cal SNIa. Values for Rv different than th e canonical Ry = 
3.1 have been found for ind ividual SNIa l|Riess et al.ll 199^ : 
Krisciunas et al. [jgOOd . l2007l ). Based on th eir entire sample . 



Reindl et all (|2005l ) find a value of 2.65. |jha et ahl (|2007l ) 



find a mean value of 2.7 for heavily obscure SN and 2.9 
for the less extincted ones. The expl anations for deviant 
Ry v alues are circumstellar material IWang fc Mukheried 
ll2007h or a different dust composition (|Riess et ahl 1 19961 : 
iKrisciunas et al.ll2006l ). Evolution in Ry may be critical for 
precision cosmology, and the SNIa themselves are also the 
best probe to characterize this evolution. 



3 EVOLUTION MODEL OF HOST GALAXY 
EXTINCTION 

Let us consider four ob servational facts: (I ) SNIa occur 
predominantly in spirals ijSullivan et al.ll2003l : iReindl et ahl 



time to rain back ( 


Howk fc Savaiid Il997l. Il999l: iHowk 


1999: Dalcanton et al. 


2004: HowkBooa Thompson et al. 


2004: iKamohuis et al 


|2007^. fill) star- formation modi- 



I2OO3I : IKrisciunas et ahl 120061 ) and (IV) there is an or- 
der of magnitude more star-formation at and beyond 
z=l (iMadau et al.lll998l:ISteidel et al.|[l999l:lGiavalisco et ahl 



I2OO4I : iThompson et al.ll2006l : iHopkind I2OO6I ) . Therefore, at 
z~l, spiral disks are likely to have an halo of ISM, recently 
processed and ejected by star-formation. For the distance 
measurement with SNIa this is especially important as their 
light will encounter more Ry < 3.1 type dust as a result (See 
Figure [T|. Hence, the average value of Ry (Ry) might be 
lower for SNIa at higher redshift. The simplest parametri- 
sation of the Ry evolution is a second order polynomial: 



Ry — az + &2 -|- 3.1, 



(1) 



with the parameters depending on the choice of Ry value 
for the different epochs. 

Of the four assertions above, the relation between star- 
formation and extinction law (HI) is the most tenuous. Ev- 
idence for a differen t extinction law in sta r-forming galaxies 
have been found bv ICalzetti et ahl (|1994| ). However, for ex- 
tended sources such as HII regions, one would need to dis- 
entangle the effects of processed dus t grains and a clumpy 
me dium (See [N atta fc Panagia) Il984l ) . 

iGordon et al. (2003) presents extinction law measure- 
ments for single stars, and hence lines-of-sight, in different 
regions in the SMC and LMC. They find evidence of dust 
processing by star-formation and lower values of Ry in the 
star-forming parts. Their average Ry values are: SMC Bar; 
Ry = 2.74, LMC SupersheU; Ry = 2.76, LMC average 
sample; Ry = 3.41. I note that the values for the regions 
with higher star-formation are lower than the average for 
the LMC. 

In this paper I explore three models: Model A assumes 
no evolution in Ry. Model B takes the LMC and SMC vah 
ues as a template for the galaxies at z=l and earlier: star- 
forming galaxies at z~l have Ry = 2.7 and gas-rich galaxies 
at z~2 and beyond have Ry »; 3.4. Model C leaves the pa- 
rameters in equation 1 free to fit the distributions of SNIa 
extinction values (Ay) in a low, intermediate and high red- 
shift sample of SNIa observations. 



4 OBSERVATIONAL TESTS WITH SNIA 

There are two suggestive observations of SNIa that support 
a general model of evolution in host galaxy extinction: (A) 
the relation between Ry for SNIa and host galaxy disk incli- 
nation and (B) the distribution of inferred extinction (Ay) 
for a sample o f high, intermediate and low redshift SNIa. 

iJha et ahl (|2007l ) use the multi-filter hghtcurves avail- 
able for a sample of 133 SNIa to independently fit the value 
of Ry as well as the light-curve peak and SNIa colour. The 
inclination and host g alaxy type are from Hyperled^ and Ry 
values from table 4 in ljha et ahl (|2007l ). Figure [2] shows the 
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Figure 2 . The relation between inclination and Ry values from 
Ijha et alj ||2007| '). for the late- type galaxies. Top panel: the num- 
ber of SNIa with Ry = 3.1 and the number with deviant values 
(dashed line). Middle panel: the percentage of SNIa with deviant 
Rv values. Bottom panel: the average value of Ry. There are very 
few SN observed in perfectly face-on galaxies. In higher inclined 
disks the average value of Ry is lower, indicating a possible effect 
of ejecta in the galaxy's halo. 



Figure 3. The cumulative histogra m of SNIa redden ing for three 
samples: the nearby SNIa sample lljha et al.l |2007|. solid hne), 
the intermediate distanc e sample from the ESSENCE survey 
(Wood -Vasev ot al. 2007', dotted line) and the high redshift sam- 
ple (Rioss et al. 2007, dashed lino). The differences in distribution 
point to evolution of the extinction in host galaxies, either in Ay, 
Ry or both. 



distribution of these SNIa as a function of disk inclination 
of the late-type host galaxies for those SNIa with Ry =3.1 
(the default) and those with deviant values {Ry 3.1). No- 
tably, the deviant values start to dominate at the higher 
inclinations. Figure [5] also shows the Ry value for each in- 
clination bin; it diminishes with increasing inclination. This 
supports our model of recent ejecta lowering the value of 
Ry . 

The second test is to compare the distributions of Ay 
of SNIa at different redshift. Figure |3] shows the distribution 
of extinction values {Ay), based on reddening and a Ry val- 
ued of 3.1, for three samples of SNIa, a local, an intermediate 
and a high-redshift on e. Local SNIa are those from the sam- 
ple of IJha et al.l (|2007l ). The intermediate redshift sample is 
the ESSENCE data from I Wood- V asev et al' ('2007') and the 
highest redshift data is from lRiess ct al. (2007). Sample sizes 
are very different: 133, 189 and 33 for the low, intermediate, 
and high redshift samples respectively. The assumption is 
that the reddening measurements are comparable between 
these samples. Each of these studies uses a 'prior' distri- 
bution for the extinction values to optimize the fit to the 
lightcurves and I assume here that the affect of this prior is 
negligible on the reported d istribution (See for a goo d dis- 
cussion of extinction priors I Wood- Vasev et al.ll2007l ). The 
distributions of extinction values in Figure[3]show that there 
is some type of evolution in Ay from one redshift sample to 
the next. SNIa searches use a "gold" and "silver" standard 
for SNIa lightcurves. To mimick these I limit the samples 
to Ay < 0.25 and Ay < 0.5 respectively in the following 
analysis. 

In Figure |4l I compare the three different models for Ry 
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Figure 4. Top panels: the three models considered for evolu- 
tion of Ry{z), (A) no Ry evolution, (B) Ry = az^ -\- hz -\- 3.1 
with Ry{z = 1) = 2.75 and Ry{z ~ 2) = 3.4, and (C) with 
Ry = az^ + fe z + 3.1 fit t o min imize the difference between the 
distribution of Ijha et al. I ll2007t) and the other two distributions 
{Rv{zmax = 0.25 to 0.3) = 3.2 to 4.1). Bottom panels: the cumu- 
lative histogram of Ay of SNIa from a nearby (solid lineljha et al.l 
[2003), intermediate (dot ted line[ Wood-Vasev ct^ .120071) ' and dis- 
tant sample (dashed line lRiess et al..2007i) for the different models 
in the top panels). 
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evolution (top panels) using the recomputed Av distribution 
for each model (bottom panels). I limit the comparison to 
the SNIa "silver" standard to maximize statistics. Figure |4] 
shows Model A (no evolution, Rv = 3.1) and two evolution 
mod els (B and C). Model B uses the Magellanic Cloud values 
from iGordon et al] (|2003h to parameterise the dependence 
of Rv{z). Model C uses the best fit of equation 1 to minimize 
the differences between the Ay distributions, assuming for 
each sample the same parent distribution, i.e., random lines 
of sight through a disk of a spiral galaxy. 

The Kolmogorov-Smirnov probability that all three Ay 
populations are from the same parent populations without 
a change in Rv (Model A, Figure |4l top left panel) is ex- 
tremely low {Pks ~ 0, Figure 131 bottom left panel). Hence, 
some type of evolution is present in the distribution of SNIa 
extinction values (Ay) or extinction law (Rv)- 

Model B is what one naively would expect on the basis 
of the narrative in §3: the star-forming galaxies have a lower 
Ry, which ris e s with gas-fraction at higher z and using the 
[Gordon et al] (|2003l ) values as a template. The top middle 
panel in Figure|3]shows the Ry values. I recomputed the val- 
ues for Ay with the model Ry{z) (Figure [4] bottom middle 
panel) . This model does not significantly improve the match 
between the three distributions of SNIa extinctions at dif- 
ferent redshifts (see Table [l|. If this model is true, there is 
substantial independent evolution in the distribution of Ay. 

Model C is a fit of the parameters in equation 1 to re- 
compute the Ay distribution and maximizing the probabil- 
ity that they are from the same parent distribution. The best 
fit values are in Table [T] and Figure [3] shows the inferred Ry 
(top right) and computed Ay distribution (bottom right). 
The best model to match the intermediate and high-redshift 
samples to the low-redshift reference peaks at Zmax ~ 0.25 
to 0.3 with a high value of Ry{z,nax)- In this case little or 
no evolution in Ay is needed. However, the values of Ry for 
each epoch are counter to expectations of §3. 

Figure [5] shows the observed average dimming of SNIa 
as a function of redshift with the maximum possible im- 
pact of Ry evolution plotted: model C with the maximum 
E(B-V) value for each SNIa allowed by the silver and gold 
sample selection. In this worst case scenario I assume that 
all SNIa are affected by the maximum reddening allowed for 
selection into the gold or silver samples, even though they 
are clearly not (Figure |3]). The worst-case scenario (pure Ry 
evolution) has the same shape as the dimming observed in 
SNIa, but the effect of Ry evolution is not strong enough to 
fully account for the observed dimming of SNIa. The exact 
origin of Ry evolution can be either changes in the imme- 
diate surroundings of SNIa or evolution in the ISM of host 
galaxies. 

Because high-redshift observations of SNIa lightcurves 
are often limited to two filters to minimize observatory time 
spent on a given object, the applied extinction correction 
is mostly the fiducial Ry = 3.1. Occasional evidence from 
higher-redshift spectra already point to very low values of 
Ry (A. Riess, private communication). Figure [5] shows that 
multi-filter and/or spectroscopic determination of Ry for 
observed SNIa are crucial to accurately determine the equa- 
tion of state of our Universe from SNIa distances. 
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Figure 5. The maximum possible difference in extinction cor- 
rection for the "silver" {Ay = 0.5 for all SNIa) and "gold" {Ay 
= 0.25 for all SNIa) samples between Ry = 3.1 and Model C. 
The points are the binned S NIa residual s, after subtracting the 
empty universe model, from iRiess et al Figure 6. Even 

if all the SNIa are affected by the maximum Ay and Model C 
holds true, evolution in Ry alone cannot explain the residuals in 
average distance moduli of SNIa. 



5 DISCUSSION AND CONCLUSIONS 

It seems likely that the distribution of Ay values for SNIa 
lightcurves will vary as a function of redshift. In this pa- 
per, I treated this as evolution of the average extinction law 
{Ry). There, however, are three possible explanations for 
the behavior of the Ay distribution, apart from a fit-prior 
artifact: 

1. Ry does not change much, but due to dust ejection 
and/or additional dust production (from the elevated star- 
formation in each galaxy), each SNIa at higher redshift is 
seen through more host galaxy dust. 

2. The host galaxies are similar to nearby spirals but SNIa 
occur preferentially in different environments at each red- 
shift. At higher redshifts, they occur in more dust-rich parts 
of spirals (e.g., the spiral arms or in interacting galaxies). 

3 Ry -and hence Ay- does evolve as a combination of 
its dependence on dust composition and dust distribution 
in the host galaxy. 

Options 1 and 2 mean that the observed difference in Ay 
distribution (Figure [4]) at different redshifts is real and the 
inferred distances from SNIa lightcurves suffer no ill effects. 
In option 3, however, the evolution in the host galaxy ex- 
tinction law could skew the observations of the acceleration 
of the expansion of the Universe and the inferred cosmology. 

I would like to stress that option 3 is not as far-fetched 
as cosmic gray dust and even the worst case scenario would 
not fully explain the SNIa dimming (Figure [5]). However, a 
realistic model would require (1) the relative positions in the 
disk of SNIa at higher redshifts to characterize the impor- 
tance of scenario 2, and (2) fits of their lightcurves which 
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Table 1. The values for two different toy models, the peak red- 
shift and peak value of Ry and the Kolmogorov-Smirnov test 
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Riess et al. (2007) 





^max 
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^max 
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gold 














(Av < 0.25) 














Model B 


1 


2.75 


0.3 
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2.75 


0.5 


Model C 


0.3 


4.1 


72 


0.25 


3.2 


97 


silver 

(Av < 0.5) 
Model B 
Model C 


1 

0.4 


2.75 
4.3 


0.2 

33 


1 

0.4 


2.75 
3.4 


2 

99 



includ e Rv as a fit parameter, similar to those of |jha et al.l 
to constrain scenario 3. Given enough statistics, the 
SNIa measurements can then be used as probes of the dis- 
tribution and severity of host galaxy extinction in disks at 
high redshift. In this paper I assumed that the distributions 
of reddening values are uniformly derived and the choice 
of prior in the lightcurve analysis did not matter. In future 
work, one would prefer the same prior in all redshift samples 
and prefer even more sufficient information to conclusively 
determine Ay and Rv for each SNIa lightcurve. 
In conclusion: 

1 There is a second plausible mechanism that modifies the 
extinction law for SNIa in higher redshift host galaxies (§3); 
host galaxy dust evolution in addition to material around 
the SNIa progenitor. 

2 The toy model of this mechanism does not match the 
Av distribution of SNIa samples at different redshifts (§4, 
Figure m Model B). 

3 Matching up Av distributions leads to a Rv evolution 
model that resembles the shape and peak of the observed 
average dimming of SNIa (§4, Figure |4l Model C). 

4 Taken as a worst-case scenario, this model still cannot 
account for the observed average dimming of SNIa with red- 
shift. 

5 It is likely that there is evolution in the distribution 
of Av values with redshift, making it impossible to deduce 
.Rv evolution from these distributions alone. As both affect 
the precision cosmology to be done with the next genera- 
tion of SNIa measurements, it is paramount that new SNIa 
observations solve for both Av and Rv separately. 

6 The bonus science of measuring both Av and Rv for 
each SNIa is a probe in the distribution, content and com- 
position of dust in their host galaxies. 
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